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Abstract 


An approach in structural optimization based 
on reliability analysis is presented with an empha- 
sis on the use of the proof-load test. In particular, 
described in detail are methods of optimizing the 
structural weight subject to a constraint on the ex- 
pected cost, which is an extended version of the 
constraint on the probability of failure. Depending 
on whether the structure is statically determinate 
or indeterminate, different methods of optimiza- 
tion are described. 

Numerical examples indicate that the expense 
of performing the proof-load test is always well 
compensated by the improvement of structural re- 
liability due to such a test. In fact, under the con- 
straint of the same expected cost, significant 
weight savings can be expected of a structure with 
proof-load-tested components, compared with the 
optimum weight of the structure consisting of com- 
ponents that are not proof- load-tested. The extent 
to which such an extra weight saving can be achieved 
depends on a parameter pertaining to the impor- 
tance of individual components relative to the cost 
of failure. 

It is shown how the significant improvement of 
statistical confidence in the reliability estimate can 
be achieved on the basis of proof-load test. The 
question of how to deal with the statistical confi- 
dence of the load distribution is also discussed at 
length. 

I. Introduction 

Thermo-mechanical properties of materials 
used for the structure of a space vehicle, such as 
fracture strength, elastic modulus, deformation 
capacity, linear thermal coefficient of expansion, 
etc. (particularly those of composite materials), 
exhibit considerable statistical variations. Fur- 
thermore, aerospace environments as well as 
loading conditions involve a number of uncertain- 
ties; for example: temperatures generated by 
aerodynamic friction, dynamic pressures, axial 
accelerations, acoustic and vibration loads, etc. 


usually tested individually or otherwise under 
simulated environmental and loading conditions 
before the vehicle is sent into the mission. Since 
such simulated tests or proof-load tests are indis- 
pensable parts of the current engineering task 
within a space program, it is extremely important 
that the effect of such tests be taken into account in 
the estimation of structural reliability and in the 
(optimum) structural design. The present study 
presents, for a given expected cost constraint, 
quantitative results of considerable weight saving 
and increased reliability by taking into considera- 
tion the proof- load test. 

From the viewpoint of reliability analysis, the 
advantage of performing the proof- load test can be 
summarized as follows. The test can improve not 
only the reliability value itself but also the statisti- 
cal confidence in such a reliability estimate. This 
is because the proof- load test eliminates structures 
with strength less than the proof load. In other 
words, the structure which passes the proof load 
test belongs to a subset, having the strength higher 
than the proof load, of the original population. 
Therefore, it is obvious that the reliability of a 
structure chosen from this subset is higher than 
that of a structure chosen from the original popu- 
lation. Furthermore, the proof-load test truncates 
the distribution function of strength at the proof 
load, hence alleviating the analytical difficulty of 
verifying the validity of a fitted distribution func- 
tion at the lower tail portion where data are usually 
non-existent. Evidently the difficulty still remains 
in the selection of a distribution function for the 
load. However, the statistical confidence in the 
reliability estimation now depends mainly on the 
accuracy of the load prediction. The question of 
how to deal with the statistical confidence of the 
load distribution is discussed in Section VI, in 
which the Baysian approach is suggested. 

The present paper: (1) develops an approach 

to an optimum design (either minimum weight 
design or minimum expected cost design) introduc- 
ing the proof load as an additional design param- 
eter, and (2) shows the practical advantage of the 
use of proof loads in terms of weight saving. 


This indicates that both strengths of a struc- 
ture and loads acting on the structure should be 
treated as random variables and that the concept of 
structural reliability should be incorporated into 
the analysis of the structure and its (optimum) 
design. In fact, some work has been done in this 
direction 1-6 at different levels of sophistication of 
reliability analysis. 

It should be observed, however, that major 
structural components of a space vehicle are 


The importance of this proposed approach in 
structural design is emphasized, not only from the 
viewpoint of optimization, but also from the reali- 
zation that it is the only rational approach in the 
face of various uncertainties and that it establishes 
a definite design procedure applicable to most 
aerospace structures. 

Although the present study places its emphasis 
on the problem of optimization of aerospace struc- 
tures, the principle involved can be applied to 
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optimization problems in other engineering 
disciplines such as the optimum design of civil 
engineering structures, naval structures, ground 
vehicles, material handling equipment and elec- 
tronic systems. In particular, the optimization can 
be highly significant, e. g. , for electronic systems 
consisting of thousands of components when the cost 
of each component is so small in comparison with 
the total cost of the system that a relatively high 
level of proof load can be applied as will be shown 
later. 

Civil engineering structures such as gantry 
towers at launching sites, buildings, bridges, etc., 
gain increasing significance for space technology 
application as, for instance, in relation to urban 
development, transportation systems and the like. 
Because of their characteristic construction pro- 
cesses, it is recognized that these structures 
usually undergo a tacit process of proof- load test 
during construction. If the structure does not fail 
during and upon completion of construction, it im- 
plies that all of its structural components and 
therefore the structure itself have sufficient 
strength to withstand at least the dead load. This 
is the information that must be taken into consid- 
eration as the lower bound of the strength distribu- 
tion for the reliability estimate of an existing 
structure. In this respect, the authors believe 
that it is not entirely impossible and in fact advisa- 
ble to devise an inexpensive method of more explicit 
proof-load test which can establish such a bound. 

Furthermore, if a structure under construc- 
tion survives a live load due to severe wind or 
earthquake acceleration, which are referred to as 
secondary live loads in many design codes but 
which are of primary importance for safety con- 
siderations of existing structures, the combined 
action of such a live load and of the dead load 
(existing at the time of occurrence of the live load) 
can be interpreted as a proof load test. The fact 
that the partially completed structure has survived 
such a proof load test should be taken into consid- 
eration in the reliability analysis since this fact 
usually makes it possible to establish a better lower 
bound of the strength of each of structural com- 
ponents (existing in the partially completed 
structure ). 

An important implication of this argument is 
that separate considerations are given to the safety 
of a structure during and after completion of its 
construction. This seems quite reasonable since 
the cost of detection possibly by means of the proof- 
load test and the cost of the replacement of that 
part of the structure which failed because of a 
member or members with insufficient strength may 
be absorbed as the construction cost, whereas any 
failure, after the structure is placed into service 
by the client, would produce much more serious 
contractual and socio-economic problems, pos- 
sibly involving human lives. 


forms are obviously possible and may be desirable 
depending on the specific problem at hand, 


EC - 


P 0i C 0i 
1 - p 0i 


+ Pf C F 


(i) 


where 


n = number of major structural com- 
ponents constituting the vehicle 

Pq. = probability of failure of a candidate 
for the i*h component under the 
stress Sq^ due to proof load. 

Pj = probability of failure of the (entire) 
structure. 


Cp - cost of failure (loss of vehicle, 
prestige, etc. ) 

Cq. = cost of proof load tests including 
cost of loss of candidate com- 
ponents that failed to pass the test. 

P 0 i/( 1 -P 0 i> “ expected number of candidates for 
the i*h component that fail under 
Sqi before the one that can sustain 
Soi is obtained. 


These quantities and hence EC are functions 
of design parameters. Furthermore, the weight W 
of the structure is also a function of the same 
parameters. Therefore, the optimization of W (or 
EC) under a constraint on EC (or W) can be per- 
formed with respect to these paramete rs. 

It is pointed out that the absolute value of Cjr 
will have no effect in the following optimization 
process. Note that if the proof- load test is not per- 
formed, then pQi = 0 in Eq. (1) and the formulation 
reduces to the minimum weight design under the 
constraint of probability of failure as discussed in 
Refs. 1-6. 

Under the further simplifying assumptions, as 
used in Refs. 1-6, that the resisting strengths R^ 

(in terms of stress such as yield stress) of indi- 
vidual components (i = 1, 2, - • v, n) are independent 
of each other as well as of the load S, the probabil- 
ity of failure p^ of the structure can be shown? to be 


Pf = 



fg(x) dx 



F r (x)f s (x) dx 


( 2 ) 


where 


H. Expected Cost of Failure and 
Optimum Design 


To present the essence of the idea, consider 
the following form of expected cost EC of the ve- 
hicle, taking only the cost of failure and of proof 
load test into account, although more elaborate 


F Ri F) 

s 

Si 

fs i < - ) 
Si 


distribution function of 
load applied to the structure 
stress acting on the i^h component 
density function of S* 

constant associated with the i** 1 
component 
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With the aid of the theoretical and experimental 
structural analysis, the stress Si acting on the ith 
component when the load S is applied to the struc- 
ture is given by 


S. 

1 


c.S 




(3) 


in which is a design parameter representing the 
size of the i th component so that its weight Wi can 
be given as Wi s? biAi [see Eq, (2 3)], ct is a con- 
stant if the structure is statically determinate 
whereas it is a function of A], • * • , A n if the 
structure is statically indeterminate, and gi(Ai) is 
a function of Ai the form of which depends on the 
nature of the ith component; for example, gi(Ai) = 

Ai = cross-sectional area for truss -like structures. 

Any method of structural analysis can be em- 
ployed to obtain Eq. (3) including the finite element 
method which is used extensively in the stress 
analysis of aerospace structures. 

The density function fSi(x) of Si can be obtained 
from the density function fs(x) of S through the 
transformation indicated in Eq. (3) and can be 
shown to be 


simplicity. It is pointed out, however, that some 
evidence exists as to its validity for statically 
indeterminate structures® as well. 


Of considerable practical importance is the case 
where the structure is subjected to a number of mu- 
tually exclusive and independent proportional loading 
systems (e.g., unusually severe wind and extremely 
strong earthquake motion) with g*(2) ( , , . ^ g^(k) 

denoting the maximum reference values of these 
loading systems within a specified period of time. 

The probability of failure of the structure is then 


where p 


(j) 


k 

given by the sum 2 ^ , WUC * C 
bility of failure ol the structure under 


P< j) 


is the proba- 
only. 


Employing in Eq, (1) the approximation indi- 
cated in Eq. (2) 


EC = 


E EC . 


( 6 ) 


i= 1 


with EC^ being the expected cost of the 


i th 


component, 


EC. = 
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1 - p 0i 


+ P fi C F 
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As was discussed in detail in Ref. 7, the fol- 
lowing points are to be noted in the derivation of 
Eq, (2): (1) the definition of structural failure is 

in accordance with the weakest link hypothesis, 
that is, the failure will take place if at least one of 
the components fails, (2) the assumption that Riare 
independent each other is a conservative one, (3) 
the approximation indicated in Eq. (2) is also of 
conservative nature, (4) the load S can be inter- 
preted as reference value of a system of propor- 
tional loading acting on the structure, and (5) if pf 
in Eq. (2) is to represent the probability of failure 
of a structure subject to a sequence of N (statisti- 
cal) loads, fs(x) should be replaced by the density 
function fg(x) of the maximum load S* in such a 
sequence; for a sequence of N “independent "loads, 
each distributed as S, the density function fsftx) is 


P 0i 

1 - p 


Oi 


( 8 ) 


with 


Y. 

i 


(«D ; 

C F 


this ratio indicates relative importance of the i** 1 
component with respect to the cost of structural 
failure and 






(5) 


where Fg(x) is the distribution function of S. 

The previous studies ^ ^ were all based on the 
weakest link hypothesis and on the assumption of 
Ri being independent of each other. The approxi- 
mation indicated in Eq. (2) was also employed in 
these studies with the exception of the work by 
Moses and Kinser , in which, essentially, the exact 
integral expression [the second member of Eq. (2 )] 
together with Eq. (5) was used in a different ana- 
lytical approach. The weakest link hypothesis 
which is probably adequate to describe the failure 
condition for statically determinate structure is 
also adopted in the present study for analytical 


p - probability of failure of the 
1 ith component 


■CO 


= 1 f r < x > f s.< x) dx 

Jo 1 1 

(9) 

The quantities EC and EC| can be expressed in 
terms of the cost of failure Cp by dividing both 
sides of Eqs. (6) and (7) by Cp 

n 

# V""' * 

EC = 2 EC. 

i=l 

(10) 

EC ? = ^ + Pfi 

(11) 
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whe re 


EC 


* 



and 



* 

which is in turn obtained from Eq. (3), where S\ 
and S* should be the measures of location of the 
same kind such as the mean value. 

It is interesting to note that upon substituting 
Eq. (14) into Eq. (4), the density function of S\ 
becomes 


Equations (10) and (11) indicate an important 
conclusion that the absolute value of the cost of 
failure Cp has no effect on the optimization pro- 
cess. In the present optimization formulation, it is 
only necessary to know or to estimate the ratio 
of the component cost Cqi to the cost of failure Cp. 

If the proof- load test is not performed, i. e. , 
if qi = 0, then EC* represents the probability of 
failure of the i*H component and EC the probability 
of failure of the entire structure. 

To perform the reliability- based optimum 
design, it is necessary to know the distribution 
function Fg^(x) of the resisting strength Ri of the 
candidate for the i*h component (henceforth re- 
ferred to as "the parent strength distribution of 
the i*h component 11 ). It is assumed that, before 
the proof- load test, this distribution and therefore 
its mean value and standard deviation are known 
with sufficient statistical confidence on the basis 
of material tests, past experience, etc. The 
question as to what is considered to be a sufficient 
confidence is crucial and will be discussed in 
detail in Part 2 of Section VI. 



and it is free from c*. 

The probability of failure poi of a candidate for 
the i^h member is given by 

»0i ■ F R.< S „ i > - F rM*) < 17 ' 

“I — l' ' 


and the probability element fR^(x) dx of the resist- 
ing strengths Ri of the proof- load- tested i*h com- 
ponent is given by 

f R (x) dx = P(x < R. < x + dx | R. > S Q .) 


Let e\ denote a design parameter indicating 
the stress level Soi of proof load to be applied to 
the candidate for the component in terms of R*; 


P(x < R. < x + dx, R. 5 S Qi ) 


0i J 


Oi 


= e.R, 
l — i 


( 12 ) 


and the central safety factor v\ be defined as 


❖ 

R. 

v= # 
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(13) 


H(x - S 0 .)f R (x) dx 

1 : F R.Ai> 

— -! 


from which it follows that 


F r _(x) 


H < X - s 0if - F R.< S 0i>] 


1 * F R. *0i 


where R* and S* are_the measures of location (such 
as the mean values R^ and S|) of the distribution of 
R^ and respectively (this notation for the mean 
value is used throughout the paper unless otherwise 
specified). 

Once the factor of safety v. i s specified, the 
candidate for the ith component should be so de- 
signed that the measure of location Si of the stress 
acting on the ith component is equal to S{ = 

Ef/ vj. This is accomplished by choosing Ai that 
"satisfies 


where 


P(E) = probability of event E 


(18) 


(19) 


P(E 1 |E 2 ) = conditional probability of Ej 
given E 2 

P(E|, E^) = probability of simultaneous 
occurrence of Ej and E^ 

H(x) = Heaviside unit step function. 


*.(A.) - 
’l' v 


c.S 

i 

❖ 
R. 
— i 


(14) Using Eq. (12) into Eqs. (18) and (19) one 

obtains 


Equation (14) is obtained from the following 

J i !_•_ 1 * T)#/ U 


relationship by replacing Sfby R*/^, 
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c.S* 
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( 21 ) 


F R.<*> 
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"(» - •**)[>«»> - F RfrE?1 


i - f 
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The Heaviside unit step function in Eqs. (20) 
and (21) truncates the original strength distribu- 
tion at the stress level of the proof load. It is for 
this reason that the distribution of R^ is referred 
to as the truncated strength distribution in the 
following. 

Since it follows from Eqs. (8) and (17) that 
qi s qi(ej) and from Eqs. (9), (16), and (21) that 
Pfi = Pfi( e i» vf), Eq. (11) can be written as 

EC* = y^ { ) + P fi (e V v.) (22) 


The optimization problem considered in the 
present study is either to minimize the structural 
weight subject to the constraint on the (relative) 
expected cost, or to minimize the (relative) ex- 
pected cost subject to the constraint on the struc- 
tural weight, both with respect to e^ and v\. Since 
the analytical technique employed here can be 
applied to either case, only the first is discussed. 
Furthermore, because the optimization analysis 
for statically determinate structures can be simpli- 
fied considerably, discussions for statically deter- 
minate and indeterminate structures are given 
separately. It should be mentioned, however, that 
the technique for statically indeterminate structures 
is a general one that can also be applied to the 
optimum design of statically determinate struc- 
tures. In fact, the optimum values in a later 
example for a statically determinate structure are 
all checked by the technique used for the statically 
indeterminate structure. 

III. Optimum Design of Statically 
Determinate Structures 


In Eq. (23), W* is the weight of the i*h 
component and bj is a known constant; for example, 
if a truss -like structure is considered, Aj is the 
cross-sectional area and b^ = 0iPi with and Pi 
being the (specified) length and density of the same 
component respectively. 

Rewriting Eq. (14), one obtains 


V. 

1 


R.g.(A.) 
—l & i' r 


c.S 

i 


(26) 


which indicates that is a function of Ai (and 
only) since cj is a constant for statically deter- 
minate structures. Hence, equivalent to Eq. (25), 


EC. = V .q.(e. ) + p-.(e., A. ) 
i n' r K fi' i i 


(27) 


where pfi(ei, A|) is used for P£i(ei, v ^(A^)) for 
simplicity. 

Now, the problem is to minimize W in^Eq. (23) 
subject to the constraint Eq. (24), with EC| given 
by either Eq. (25) or Eq. (27). 

Since Eq. (23) is linear and there is only one 
constraint equation, Eq. (24) is always active; 
i. e. , the equality sign of Eq. (24) holds at 
optimum. 

The present problem can now easily be formu- 
lated with the aid of the variational principle. 
Although the method does not generate the solution 
explicitly, it indicates that, for a minimum weight 
design, the following relationships must be satis- 
fied, under the assumption that EC* is small com- 
pared with unity so that at optimum, the variation 
of Wj/W is small in comparison with that of ECj 
(see Appendix). 


The optimization problem in this case can be 
stated as follows. 

Minimize the weight of the form 


OEC. 

l 

ae. ' 

1 


= 0 


(i = 1, 2, 


, n) 


(28) 


W = 


a a 

E 'V" ip % 

w i = Ev, 


i=l i=l 


(2 3) 


EC* W. 
1 1 


ec; 


= 1 r (i = 1, 2, 


, n) 


(29) 


subject to 


with 


EC = 2^ EC i 5 EC a 
i=l 


EC. = V i q i (e.) + p fi (e i , V .) 


Equations (28) state that, for an optimum 
structural weight, the stress level of the proof 
load to be applied to individual components should 
also be optimum in the sense that the relative cost 
EC? of individual components are minimized at 
(24) that stress level. As an example, under the 

assumption that R^ and Sj are normally distributed 
with coefficients of variation of 20% for Sj and 5% 
for R*, the dependence EC* on e^ using a specific 
value of v i (=1. 6) is plotted in Fig. 1 with Yi as a 
parameter. The locus of those points at which 
ECf assumes minimum values (Curve I) plays an 
(25) important role in the following optimization proc- 
ess. Those values of the proof load ei that produce 
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a Function of Proof Stress Level ei; Load Distri- 
bution is Normal with Coefficient of Variation 
0. 20 and Parent Strength Distribution is 
Normal with Coefficient of Variation 0. 05 


minimum EC* associated with given vj and Yi are 
denoted by e*. This implies^that Curve I indicates 
the relationship between EC* and e* given and 
Vi- 

Equations (29) state than an optimum weight is 
realized when the total weight is allocated to indi- 
vidual components in proportion to their expected 
costs. This fact was shown to be valid also for 
optimization without proof-load test^ in which case 
the total weight was to be allocated proportionately 
to the probabilities of failure rather than to the 
expected costs. 

Usefulness of Eqs. (28) and (29) lies in the 
fact that these can be used to develop an iterative 
procedure consisting of the following steps in 
arriving at a minimum weight design: 

(1) Construct a diagram in which the 

ECf ~ ef relationship is given for various 
values of v\ and Yi (Fig. 2). This is an 
extended version of Fig. 1. 

(2) Try EC* = EC*/n as a first estimate of 
EC*. 

(3) Read, from Fig. 2, ef and correspond- 
ing to the latest estimate of ECf and to 
the specified value of Yi* Aj is then 
calculated from Eq, (14). 

(4) Compute EC* from Eqs. (29) using those 
values of A i just obtained in (3). 

(5) Go to (3) with the value of EC* estimated 
in (4) and repeat the procedure. 



Fig. 2. Relative Expected Component Cost EC^ as 
a Function of Optimum Proof Stress Level ef; Load 
Distribution is Normal with Coefficient of Vari- 
ation 0.20 and Parent Strength Distribution is 
Normal with Coefficient of Variation 0. 05 


The rate in which this process converges to 
stable values of W j and ei is extremely rapid since 
the component weight Wi is insensitive to the vari- 
ation of EC*. In fact, experience shows that two 
cycles of iteration are sufficient to obtain the 
optimum design. 

It is to be noted that in a numerical example given 
later, the information contained in Fig. 2 is stored 
in the computer memory and the third step of the 
above procedure is accomplished by means of an 
interpolation subroutine. 

IV. Optimum Design of Statically 
~~ Indeterminate Structures 

For statically indeterminate structures, the 
central safety factor s v{(i = 1 , 2, - * • , n) cannot be 
chosen arbitrarily since they should satisfy con- 
tinuity equations. This situation can readily be 
demonstrated by a statically indeterminate truss 
subjected to a system of proportional loading with^ 
a reference value S. The applied stress Si (and Sj) 
acting on individual members can be determined 
only after the cross-sectional area Ai of the mem- 
bers are specified. This implies that vi are func- 
tions of Ai, A'2, • • • , A n . In fact, ci in Eq. (26) 
are, in general, functions of A i, Az, •••, A n . 

Hence, 


v. = v (A . , A , • * * , A ) i = 1, 2, * * • , n (30) 

li i c. n 


It is important to note that the inverse of Eq. 
(30) does not exist. In other words, one cannot 
express Ai as a function of vj, » ■ * * » anc ^ 
therefore vj(i = 1, 2, • • • , n) cannot be used as inde- 
pendent design variables as in the preceding section 
The optimization problem should therefore be stated 
as follows: 

Minimize 


n 

w-SVi (31) 

i=l 
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subject to 


ec * = 22 ec * - ec I 

i=l 


(32) 


with 


EC i = + Pfi< e i« ^ 


(33) 


or 


EC i = Vjq^e.) + p fi (e i ,A 1 ,A 2 ,---.A n ) 


(34) 


The basic difference between the present case 
and that of statically determinate structures is 
quite clear; the probability of failure pfi now de- 
pends not only on A* and e i but also on A \ , 

Az, ••v, Aj_|, Ai+i, • • •* A n . Unfortunately this 
fact makes it impossible to solve the present opti- 
mization ("design") problem in a simple iterative 
approach efficiently employed for statically deter- 
minate structures in the preceding section. It is 
pointed out in passing, however, that such a dif- 
ference between statically determinate and indeter- 
minate structures hardly makes it necessary to 
treat these structures differently in the reliability 
"analysis" where the essential problem is to esti- 
mate the probability of failure of "designed struc- 
tures" as long as the weakest link hypothesis is 
assumed regardless of structural determinancy. 

Nevertheless, the variational principle, as 
described in the Appendix for statically deter- 
minate structures, can be applied to the problem 
of indeterminate structures at least formally to 
obtain the following conditions for an optimum 




3EC. 

i 


d e . 

l 


0 


(i * 1, 2, * • • , n) (35) 


Hence, the problem is now to minimize W in 
Eq. (31) by a proper choice of A{ and e£ subject to 
a constraint on the expected cost Eqs. (34) and sat- 
isfying Eqs. (35). 

It is believed that the optimization technique 
most appropriate to the present problem is a 
gradient move method as briefly described below. 

In this method, a design point Bj is first 
chosen arbitrarily in the acceptable domain defined 
by EC* < EC* of the n-dimensional design space of 
Al, AZf *••.# A n . 

Note that once Bj is chosen, Aj f A2, •••* A n 
are given and v% can be computed from Eq. (26). 
With these values of vi and the specified values of 
Vj, ECf and e* can be read from Fig. 2. This 
makes it possible to check if B] is in the acceptable 
domain. 


The design is then modified by moving normal 
to the weight contour by a specified step from point 
B] to a new design point B 2 with a lighter weight. 
This process is repeated until the constraint EC*= 
EC* is reached at point Bq. 

Let U and V be respectively the gradients of 
the relative expected cost EC* [with e{ replaced by 
optimum values ef satisfying Eq. (35)] and the 
weight W at point Bq (and hence normal to the con- 
straint EC* = EC^ and the weight contour); 




k=l 


_ • V«* 

i k X k 

(36) 


v-w,£|w lk .£ 


k=l 


Vk 


(37) 


k = 1 


with i k being the unit base vector in the positive 
direction of A^ axis. 

Let Q be a vector such that 


(U- 

Q, “ b o S 0 

(38) 

(V- 


(39) 


The direction of Q defines the so-called usable 
feasible direction. 9 A systematic scheme for find- 
ing Q has been proposed by Zoudendijk^ and used 
in the present study. 

Since EC^ are functions of and , the 
partial derivatives 


dEC* 

1 




at B- 


(i = 1,2, ••*, n) 


can be obtained from Fig. 2 by interpolation, 
whereas the partial derivatives 


BA k at B q 


(i = 1,2, • • • ,n, k = 1, 2, • • • , n) 


from Eq. (26) with the aid of the finite difference 
technique. 


The design point is now moved from Bq along 
Q in a specified step away from the constraint 
EC = EC:£ into the acceptable domain with a reduc- 
tion in the weight. The modification of the design 
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proceeds along Q until the design point reaches the 
constraint again. Then, another usable feasible 
direction is found and the process is repeated until 
the design point B* is reached on the constraint at 
which the Kuhn-Tucker optimality condition*? is 
satisfied [Q cannot be found at B* so as to satisfy 
Eqs. (38) and (39) simultaneously with at least one 
of them being purely an inequality]. This point B* 
is an optimum design point corresponding to a local 
minimum for the weight W. The global minimum 
can usually be found as the least of the local minima 
obtained by beginning with a number of starting 
design points. 

It should be mentioned that disregarding Eqs. 
(35), a straightforward application of the gradient 
move technique can be made by taking ei(i = 1, 

2, ■•••-, n) as independent design variables in the 
design space. Hence, a design point B i is first 
chosen arbitrarily in the acceptable domain defined 
by EC* < EC| of the 2n- dimensional design space 
of At, A 2 , ** *» A n , 6 j, ®2* * ' * > fin- Then, the 
procedure just described can be employed to obtain 
a local minimum. It is believed that the compu- 
tational work involved in such an approach will be 
much more than the approach taking advantage of 
Eqs. (35). It is emphasized that Eqs. (35) provide 
not only the computational advantage but also a 
physical significance of the optimum test level of 
components, as discussed previously in Section III. 

V. Numerical Examples 

1. For the purpose of comparison, the same 
numerical example as in Ref. 6 is considered in 
which the minimum weight design is performed of 
a statically determinate, truss -like structure con- 
sisting of ten components and subjected to a system 
of proportional loading. Since the failure condi- 
tion is assumed to be that of yielding, the resisting 
strength Rj = try = yield stress. The assumption of 
Vi being a constant for all components is used for 
simplicity without loss of generality. 

A constraint imposed on EC in this example is 

EC £ 10" 3 C„ or EC* < 10" 3 (40) 

r 

Note that, without the proof load test, this 
formulation reduces to the minimum weight design 
under the constraint of probability of failure 
p£ S 10“ 3 as discussed in Ref. 6. 

The weight W can be written in the form of 
Eq. (2 3) with n - 10 and 

pL.c .H v. 

W. = pL.A. = — kJ- — - (41) 

1 11 <T 

y 

where 

p = density of material = 0. 283 lb/in. ^ 

L.^ = length of the ith component = 60 in. 

S = mean value of S = 60 x 10 3 psi 

7r = mean value of yield stress 

Y o-y = 40 x 10 3 psi 

c,, = 0. 1 x i(i = 1, 2, • • • , 10) 


It is assumed that the distribution functions of 
S and (Ty are both normal with the coefficient of 
variation 0. 20 and 0. 05 respectively. Then, a choice 
of a set of values for v{ and Soi determines the 
distribution function of Rj (that is, the truncated 
distribution of <ry). This in turn makes it possible 
to evaluate pfi [Eq. (9)]. Therefore, in this 
example, the independent design parameters are 
v£ and ej or Soi - e i&i “ e i <r y^ 

The minimum weight design subject to the con- 
straint of Eq. (24) can now easily be achieved by 
the iterative method described in Section HI. 

The result is shown in Table 1 where the con- 
straint of pf £ 10*3 i s used for both conventional 
design and standard optimum design (without proof 
load test). Since all these designs are associated 
with the expected cost of 10-3 Cp, Table 1 indi- 
cates the fact that, by performing the proof- load 
test, not only the reliability of the structure in- 
creases, but also considerable weight saving is 
achieved. 

It is further observed from Table 1 that the 
extent of weight saving and reliability increase 
depend essentially on the value of Vf, which is the 
ratio of the cost of the ith component with respect 
to the cost of failure. For smaller values of Vf, 
larger proof loads can be applied with the same 
constraint of the expected cost, thus yielding lighte r 
structural weight (this can easily be realized from 
Fig. 2) and smaller probability of failure. For 
example, the structural weight W and the prob- 
ability of failure pf associated with the optimum 
structural design without proof load test are 25 3. 2 
lb and 10“ 3 respectively, while those associated 
with the optimum structural design with proof load 
test are 221 lb and 0. 613 x 10-3 for V| = 10-6, and 
243. 9 lb and 0. 625 x 10-3 for Vj = 10-4. 

Therefore, as a result of the proof- load test, 
one concludes that, for the optimum design, higher 
benefit can be obtained for smaller values of Vi. It 
is due to this conclusion that the optimum design 
with proof load test can be highly significant for 
electronic systems consisting of thousands of com- 
ponents when the cost of each component is very 
small compared to the total cost of the system. 

The study in Ref. 6 showed that more weight 
saving can be expected if the more accurate expres- 
sion rather than the approximation is used for 
evaluating the probability of failure [Eq. (2)], 
although then Eqs. (29) will no longer be valid 
and hence a more elaborate computational scheme 
has to be applied. 

2. A (statically indeterminate) three -member 
truss is designed for a minimum weight to resist 
a set of proportional loading as shown in Fig. 3. The 
mean value of yield stress for each member is 
40 x 103 psi and that of S is 100 x 103 psi. The 
constraint on the cost EC is 5 x 10-4 Cp. Both tr y 
and S are normal with coefficients of variation 
0. 20 and 0. 05 respectively. 

The weight function is 

^ = (2) 1/2 A x + A 2 + (2) 1/2 A 3 (42) 

where W = total weight, p = density of material, 

L. = length of member 2 and Af = area of the ith 
member. 
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Table 1. Ten- Member Structure (EC^ = 10'**) 




Conventional 
design (equal 
safety factor) 

■ 

Standard 

optimum 

design 

Current optimum design (with proof- load test) 

Member i 

Mean load, 
S 

Y. = 

l 

io - 6 

Y = 

10' 5 

■ 

IO' 4 

V. = 10" 3 
1 




m 

m 

c* 

i 

l 

* 
€ , 

1 

Ai 2 

in.t 

* 
€ . 
1 

I5S 


1 

0. 1 

0.274 

0. 287 

0. 257 

1. 064 

0.271 

0. 998 

0.283 

0. 920 

0. 287 

0.84 

2 

0.2 

0. 547 

0. 562 

0. 498* 

1.079 

0. 524 

1. 019 

0. 550 

0. 941 

0. 561 

0. 86 

3 

0. 3 

0, 817 

0. 833 

0. 734 

1. 079 

0.771 

1. 031 

0.812 

0. 954 

0. 831 

0. 87 

4 

0.4 

1, 09 

1. 101 

0. 966 

1. 087 

1.013 

1.D39 

1. 068 

0. 963 

1. 097 

0.88 

5 

0. 5 

1. 37 

1. 367 

1. 196 

1. 096 

1.252 

1. 044 

1. 322 

0. 971 

1. 361 

0. 88 

6 

0.6 

1.64 

1. 630 

1.424 

1. 100 

1.490 

1. 049 

1. 573 

0. 977 

1. 622 

0. 89 

7 

0.7 

1.92 

1. 893 

1.65 

I. 103 

1.725 

1. 053 

1. 821 

0. 982 

1. 882 

0. 89 

8 

0. 8 

2. 19 

2. 153 

1. 875 

1. 105 

1.959 

1. 056 

2. 068 

0. 986 

2. 140 

0. 90 

9 

0.9 

2.46 

2. 413 

2. 098 

1. 107 

2. 192 

1. 058 

2. 313 

0.990 

2. 397 

0.90 

10 

1. 0 

2.74 

2. 672 

2. 320 

1. 109 

2.423 

1. 061 

2. 556 

0. 993 

2. 653 

0. 91 


Structural 
weight 
(W) lb 

255.6 

253.2 

221. 0 

231. 3 

243. 9 

251 

. 8 


Probability 
of failure 

(PC) 

1. 0 x 10- 3 

1.0 x 10‘ 3 

0.613 

X 10-3 

0.6158 

x 10- 3 

0.625x IO- 3 

0.839 x IO- 3 



The gradient move method described in the 
preceding section is employed to find an optimum 
design. It is found that, under the loading condi- 
tion described in Fig. 3, the minimum weight 
design is the one for which the area of member 1 
is zero. Or, in this particular case, the minimum 
weight design is a statically determinate structure. 
For the purpose of comparison, the result is listed 
in Table 2 for different values of V- (Y^ are as- 
sumed equal for i = 1 , 2, 3) as well as for zero 
proof load (standard optimum design). Again, 
considerable amount of weight saving is accom- 
plished for smaller values of Y^. 


Although the possibility of buckling is not con- 
sidered here, it can be treated without any 
difficulty. 7 

3. A spherical shell of constant thickness 
fixed around its edge and subjected to a uniformly 
distributed load S is to be designed for a minimum 
weight (Fig. 4). The mean values of the yield 
stress o*y and the applied load S are respectively 
<Ty = 45 x 1()3 psi (for both tension and compression) 
and S = 0.6 x 103 psi. The constraint on the ex- 
pected cost EC is EC a = 10-5 Cf* Both cry and S 
are normally distributed with coefficient of varia- 
tion 0. 05 and 0. 2 respectively. The maximum 
stress <r max due to load S is the meridional stress 
at the fixed edge and approximately equal to i0 

o- = -Sf(a, o,h) (4 3a) 


with 

f(a, a, h) = |iLJLi£LfijTo, 7 5 - 0. 038 sin 2 * 


for 

a -gjg— - < 3 (43b) 

h 
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Table 2. Three-Member Truss (EC* n 0. 5 x 10“ 3) 


Member 

Standard 

Optimum 

Design 

11 

10- 6 

v = 

10" 5 

Y. = 
1 

10~ 4 

= 

10’ 3 


in.^ 

Ai 2 

* 

e. 

1 

H 

s 

e. 

1 

■ 

* 

e. 

1 

■ 

❖ 
e , 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

4. 568 


1. 119 

4. 102 

1. 07 5 

4. 320 

1, Oil 

4. 515 


3 

1. 579 

1 

1. 107 

1. 433 

1. 059 

1. 512 

0. 991 

I. 568 

mm 

w 

pL 

6. 802 

5. 

88 

6. 129 

6. 

459 

6.732 

Pf 

5 x 10" 4 

3 x 

10- 4 

2. 8 x 

0 

1 

4* 

3. 2 x 

h— 

O 

1 

4* 

3. 8 x 

10’ 4 


S 



With the aid of Eqs. (4 3), the mean applied 
stress <r max can be defined as the maximum stress 
produced by the mean applied load S; o- max = S 
f(a,j>, hj. The central safety factor v is then 
v = ^y/^rnax jn-d the stress level Sq due to proof 
load is Sq = ec r y . 

Since this is a one-component structure, the 
optimum design can be achieved without using the 
gradient move method. The procedure is as follows: 

1. Construct a diagram where the EC*- e* - vj 
relationship is given for various values of Yi 
(Fig, 2). 

Z. Read e* and v from Fig. Z for specified 
constraint EC| and given value of Y. 

3. With the safety factor v just evaluated, the 
thickness h of the shell is computed using the fol- 
lowing expression obtained from Eqs. (43). 


f(a, oyh) = (44) 

vS 


and 


f(a,«,h) = 1.2^ 


for 


. 2 

asm a 


< 1Z 


where a is the shell radius. 


The results are listed in Table 3 for various 
values of Y including the case of standard optimum 
design (e* = 0) for the purpose of comparison. 

4, A spherical shell with variable thickness 
and subjected to a uniformly distributed load S is 
designed for a minimum weight (Fig. 5). The 
thickness at the top of the shell is h 1 while it is h 2 
at the clamped edge. The thickness varies linearly 
with respect to the angle o. No attempt is made 
here for the determination of the optimum shape 
of the shell (for this aspect, the reader is referred 
to Ref. 11). The following values are used for 
numerical computation; S = 380 psi, cry = 40 x 10 3 
psi, E = 30 x 10& psi, G = 12 x 10& psi, p = 0. 2 38 
lb/in. ^ Both cry and S are assumed to be normally 
distributed with coefficient of variation 0. Z and 
0, 05 respectively. 
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Table 3. Spherical Shell with Constant Thickness 
(EC* = 1. 0 x 10-5) 



Standard 

optimum 

design 


r — — — ? — 

y = io " 5 


V = 10 3 

h, in. 

0. 315 

0. 28 

0. 297 

0. 31 

0. 314 

e* 

0 

1. 05 

0. 983 

0. 895 

0. 824 

V 

1.971 

1.75 

1. 854 

1. 941 

1.964 

Pf 

1 x 10' 5 

5. 9 x 10' 6 

5. 6 x 10“ 6 1 

6.7 x 10' 6 

7.2 x 10' 6 



Fig. 5. Spherical Shell of Variable Thickness 


The weight function is 



(45) 


where a is the shell radius. The finite element 
method is employed for stress analysis of the shell 
while the gradient move method is used to find the 
optimum design. The maximum stress occurs 
either at the clamped edge (meridional stress) or 
at the top of the shell (tangential stress) depending 
on the magnitude of h\ and h^. 

It is noted that a number of starting design 
points were tried all resulting in the same optimum 
design. The results are listed in Table 4 for dif- 
ferent values of Yj and EC*. For an easy refer- 
ence, optimum values of thickness (hi and h?) and 
the weight (W) are plotted as functions of EC^ in 
Figs. 6 and 7, respectively. 

VI. Discussion 


1. The result of the preceding examples indi- 
cates, as expected, that the level of proof load to 


be applied to individual components is lower for 
more important components with a larger value of 
and higher for less important components with a 
smaller value of Yj, reflecting a simple fact that 
the more expensive the component is, the less one 
can afford a possibility to lose it by the proof load 
test. 


For instance, in the truss considered in 
Example 2, under the same constraint of EC* = 

5 x 10-4, the optimum levels of proof load to be 
applied to member 3 are e* = 0. 907 for Y = 10-3 
and e* = 1. 107 for Y = 10-6. Similarly, the opti- 
mum values of e^ are 0. 824 for Y = 10-3 and 1. 05 
for Y = 10-6 i n the spherical shell of Example 3. 


2. The statement is made in the Introduction 
that the proof-load test can improve the statistical 
confidence in the reliability estimate because the 
test truncates the distribution function of the 
strength, hence alleviating, if not completely re- 
moving, the difficulty of justifying the use of a 
fitted distribution function at the lower tail portion 
where data are usually non-existent. 


The validity of this statement evidently rests 
on whether the truncated strength distribution can 
really be established with a significantly improved 
confidence on a sample of practical size. In gen- 
eral, this can be achieved if the magnitude of the 
proof load is reasonable in the sense that it is 
equal to a strength value within a central portion of 
the parent strength distribution, since then the sig- 
nificant part of the truncated distribution in con- 
nection with the evaluation of probability of failure 
(strength values larger than but close to the point 
of truncation or the proof load) involves neither 
extreme lower nor extreme upper tail of the parent 
strength distribution. Note that this is equivalent 
to a statement that p q| should not be too small com- 
pared with or too close to unity. For example, 
when considering a coefficient of variation of parent 
strength distribution of 0. 05, values of e^ in the 
±2<r range (between 1. 0 - 2 x 0. 05 = 0. 90 and 
1.0 + 2 x 0. 05 - 1. 10) may be regarded as reason- 
able. Hence, the optimum levels of proof load test 
obtained for member 3 in Example 2 are reasonable. 
However, the optimum level for Y = 10“ 3 in Exam- 
ple 3 is not reasonable because the proof load is so 
small that a sample of unreasonably large size 
would still be required to establish the strength 
distribution although truncated by the proof load. 
Similar situation exists when the level of proof load 
is unreasonably high. 
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Table 4. Optimum Design (hi,h2 in in., W in lb) 


Y - 10- 6 

EC* 

0. 3 x 10' 4 

0. 3 x 10' 3 

0. 3 x 10' 2 

0. 3 x 10' 1 

h l 

0.755 

0.698 

0. 626 

0. 547 

- 2 

1.694 

1.532 

1. 372 

1. 187 

V 

1.678 

1. 511 

1. 349 

1. 171 

e* 

1. 07 3 

1. 112 

1. 149 

1. 174 

w 

5798 

5266 

4716 

4078 

Pf 

►— » 
O 

o 

X 

o 

1 

U1 

2.22 x 10' 4 

1. 62 x 10' 3 

2.64 x 10' 2 

Y = 10' 5 

EC* 

0. 3 x 10" 4 

0. 3 x 10' 3 

0. 3 x 10' 2 

0. 3 x 10" 1 

h l 

0. 804 

0.7136 

0.647 

0. 561 

h 2 

1.791 

1. 596 

1.415 

1.221 

V 

1.768 

1. 573 

1. 392 

1,203 

e* 

1. Oil 

1. 072 

1. Ill 

1. 145 

w 

6139 

5465 

4867 

4205 

Pf 

1. 59 x 10' 5 

1 

1. 76 x 10' 4 

2. 26 x 10“ ^ 

2. 46 x 10“ 2 

■ 

II 

o 

i 

4^ 

EC* 

0. 3 x 10' 4 

0. 3 x 10' 3 

0. 3 x 10' 2 

0. 3 X 10' 1 

h l 

0. 85 

0.7 52 

0. 665 

0. 581 

*2 

1.88 

1.672 

1. 474 

1.264 

V 

1. 855 

1.657 

1.45 

1.242 

e* 

0.932 

1. 009 

1. 068 

1. 108 

w 

6451 

57 32 

5058 

4351 

Pf 

1. 99 x 10' 5 

1. 67 x 10' 4 

1. 96 x 10' 3 

2. 38 x 10' 2 

Y = I0' 3 

EC* 

0. 3 x 10" 4 

0. 3 x 10' 3 

0. 3 x 10' 2 

0. 3 x 10' 1 

h l 

0. 859 

0.796 

0.700 

0. 603 

- 2 

1. 9 12 

1.7 53 

1. 548 

1. 316 

V 

1. 885 

1.7 30 

1. 527 

1. 296 

e* 

0. 857 

0. 927 

1. 003 

1. 061 

w 

6554 

6021 

5314 

4529 

Pf 

2. 76 x 10' 5 

2. 23 x 10' 5 

1. 90 x 10' 3 

2. 19 x 10' 2 

Standard Optimum Design 

EC* = Pf 

0. 3 x 10' 4 

0. 3 x 10' 3 

0. 3 x 10“ 2 

0. 3 x 10' 1 

h l 

0. 860 

0.80 

0.728 

0. 650 

h 2 

1.916 

1.771 

1. 616 

1.42 1 

V 

1. 888 

1.749 

1. 592 

1.399 

e* 

0 

0 

0 

0 

w 

6566 

6078 

5541 

4886 
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AND h„, in. 



Fig. 6. Optimum Thickness hj and of 
Spherical Shell of Variable Thickness vs 
Expected Cost Constraint EC* 



Fig. 7. Optimum Weight W of Spherical Shell 
of Variable Thickness vs Expected Cost 
Constraint EC* 


This observation makes it necessary to empha- 
size the tacit assumption employed in the preced- 
ing examples that the cost CEi of the test for 
establishing the truncated strength distribution of 
the i*h component (this test should not be confused 
with the proof-load test) is independent of the point 
of truncation (ej or Soi) yet to be determined. Such 
an assumption, employed in the present study for 
first approximation, is believed to be valid when 
the statistical properties of the strength of the 
material used is well known on the accumulated 
data or when the optimum levels of the proof load 
are found to be within the ±2<r range of the parent 
strength distribution as a result of the analysis. 
Evidently, the cost Cjri(ei) of the test for estab- 
lishing the truncated strength distribution with a 


reasonable confidence is expected to increase 
rapidly as the optimum location of the truncation 
increases or decreases beyond the ±Z(T range, par- 
ticularly when the material to be used is new to the 
engineers. 

It is strongly emphasized that the assumption 
°f C Ei being independent of ei is not essential in 
the present formulation and analysis. In fact, the 
effect of the cost Cjri can easily be included in the 
formulation by modifying Eq. (I) into the form 


EC ’ Z * C Ei , 'i ) ] * »( C F 


provided, of course, that the functional form of 
^Ei( e i) is reasonably well known. 

The formulation based on Eq. (46) is much 
more likely to produce optimum values for ej 
within a central portion of the parent strength dis- 
tribution because other values of e{ incur extremely 
large cost in CEi( e i) anc ^ therefore are unlikely to 
become optimum values. 

This fact is of utmost importance. Indeed, it 
is because of this that the optimum design based on 
reliability and proof-load test is expected to be 
insensitive to the tail portions of the parent strength 
distribution. 

Even in the formulation without involving the 
cost Cjri, most of the cases demonstrated in 
Examples 1-4 produced reasonable optimal levels 
of proof load implying that most of these optimum 
solutions were insensitive to the analytical form of 
the parent strength distribution. 

This is the main reason why the normal distri- 
bution was assumed for the parent strength distri- 
bution without apparent justification at the outset. 
Another reason is that all of the previous work also 
assumed the normal distribution for the parent 
strength distribution and therefore the same had to 
be assumed in the present study for possible 
comparison. 

It is recommended, however, that the Pearson 
distribution family 12 be used if the first four 
moments of the distribution are known reasonably 
well. 

3. In the present approach, the strength within 
an individual member is assumed to be invariant 
(though statistical) for simplicity. Such an assump- 
tion is, however, subject to a critical observation 
that in reality the strength usually varies statis- 
tically from point to point even within the same 
member. It should be recognized, therefore, that 
the failure of the member actually occurs once the 
resisting strength at any location within the member 
is exceeded by the load (in terms of stres s , strain, 
or displacement depending on the definition of 
failure) acting at that point. For example, when 
the load is uniformly distributed, as in a truss 
member, the resisting strength Ri of the i*h mem- 
ber can be defined as the strength of the weakest 
cross-section of the member and, therefore, may 
be distributed according to one of the asymptotic 
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distribution functions of the smallest values. This 
implies a possible necessity of taking the statis- 
tical size effect into consideration for a more 
rigorous analysis. 

It is pointed out that a similar but possibly 
more complicated situation arises if, for example, 
elastic moduli of the material are more realis- 
tically treated as statistical functions of spacial 
coordinates. In such cases, the coefficients of 
the constitutive equations become statistical func- 
tions of spacial coordinates, hence making the 
stress analysis untractable, at least at present 
time, even under the assumption that the statis- 
tical functions are homogeneous. 

4. The normal distribution is used for the 
load as an example simply because, as for the 
strength distribution, the same was used in the 
previous work. The use of the normal distribu- 
tion is, however, not essential for the development 
of the present analysis. Therefore, other distri- 
bution functions should be used if there is any 
reason to believe empirically or theoretically that 
they represent the statistical load better than the 
normal distribution. 

The following discussion proceeds under the 
assumption that the load is normally distributed. 
However, the discussion applies in principle to the 
case where the load is distributed otherwise. 

It is only in rare occasions that even a rea- 
sonable amount of data exist for the environmental 
condition of a specific space mission. For the 
prediction of the load therefore, engineers usually 
must depend on incomplete knowledge and past 
experience if such has been accumulated on similar 
missions. Under these circumstances, the reli- 
ability of the two parameters (the mean value ps 
and the standard deviation <rg) of the normal distri- 
bution assumed for the load suffers from a con- 
siderable lack of statistical confidence. 

One possible way to cope with this kind of 
situation seems to be a sensible use of the Baysian 
approach in which these parameters are treated as 
if they were random variables with a joint density 
function 4>(x, y) (x for ps a *id y for ^S) constructed 
so as to reflect both the past experience and the 
accuracy of the load prediction. If a set of ob- 
served data directly related to the load are some- 
how available, such information should be used to 
modify the density function 4>(x, y), following the 
Bayes theorem, by multiplying it by the likelihood 
of the observed data. 

It is now clear that the minimum weight W* 
computed in the preceding section is a conditional 
one under a given set of ps and <rg; W* = W*(pg,ag). 
In other words, depending on the values of the 
parameters pg and (rg, the minimum weight 
assumes different values. It is important to real- 
ize that W*(pg,crg) can be interpreted as n the best 
design if the mean and the standard deviation of 
the load distribution are truly ps and <r g. M 

The question immediately arises then which 
design one should choose in the face of uncertainty 
involved in the mean value and the standard devia- 
tion. A possible answer to this question is as 
follows. First, construct a loss function that 
represents, in some general sense, the loss L 
caused by the choice of a specific minimum weight 


design Wq in place of W (pS»<^s) that should be 
chosen if the mean and the standard deviation were 
known to be ps and erg respectively. 

The analytical form of the loss function would 
probably depend on the managerial as well as engi.- 
neering judgment except for the fact that it is a 
function of | Wo - W*(pg,<Tg)|. 

Once the form of the loss function is con- 
structed, then "a best design in the face of 
uncertainty on the mean and the standard deviation” 
is chosen as the design that minimizes the expected 
loss EL or the expected value of L with respect to 
pg and org. In other words, satisfies 



with D being the two dimensional domain in which 
x and y are defined. 

For example, if a simple quadratic loss is 
assumed for L as a first approximation, 

L = [ w 0 - * V ' s >] 2 (49) 

without discriminating an under-weight design 
against an over-weight design, then it follows 
immediately that 

wj = E [wV. s )] (50) 

To compute the^ expected value in Eq. (50), the 
minimum weight W*(ps,crg) has to be evaluated at a 
reasonable number of sets of values of ps and erg. 
This implies that the optimization procedures 
described in the preceding sections must be re- 
peated the same number of times. The cost of per- 
forming such computation, however, may or may_ 
not be justified. For practical purposes, W*(pg,crS) 
may be a good approximation for E[W*(ps,o'S)] 
although it is not quite clear at this time how accu- 
rate this approximation is. 

In this context, the result of the preceding __ 
numerical examples can be considered as W*(pSiO*s) 
if the numerical values used for the mean and the 
standard deviation of the load distribution are inter- 
preted as the best estimate pg (denoted by S in the 
preceding sections ) and <rg of pg and o-g respectively. 

It is pointed out, that such an approximation 
is in general reasonable in view of the fact that the 


14 



analytical form of the loss function itself is a 
product of subjectively inclinded engineering, 
economical and managerial judgment. 


An alternative approach of choosing design is 
that of the Baysian confidence as described below. 


Consider, for example, 
values pgp and <rgp such that 


a set oi 


pai aiuc ic i 


described in detail. Different methods of optimi- 
zation, depending on whether the structure is 
statically determinate or indeterminate, are de- 
scribed. The formulation of optimization problems 
using the constraint on the expected cost as defined 
in the present study is more general than using the 
constraint on the probability of failure as employed 
in the existing literature in the sense that the 
former reduces to the latter if no proof- load test 
is performed. 


P(p s < M-gp# o-g < Ogp) - p (51) 


Furthermore, consider a design with the minimum 
weight W*(h>SP» o'SP) assuming that psp and ar$p are 
the mean value and the standard deviation of the 
load distribution. Then, the probability that the 
best minimum weight ^S) associated with 

(unknown) true parameter values ps> will be 
less than W*(fiS(3, ^SPK is P* It is assumed in the 
above statement that the smaller the mean value 
and the standard deviation are, the smaller the 
resulting minimum weight is. A design associated 
with jjisp and ogp can be considered reasonable if 
the confidence coefficient is small. It is to be 
noted that in this approach, the design becomes 
more conservative as a smaller value of p is 
specified. 

5. At the last stage of this investigation, a 
paper by Barnett and Hermann 13 came to the 
attention of the present authors. It is acknowledged 
that the paper recognizes the practical significance 
and importance of proof testing from the viewpoint 
of reliability analysis and suggests a method of 
component optimization on the basis of theWeibull 
strength distribution [see part 3 of this section in 
this respect]. 

The present paper, however, explicitly 
describes methods of reliability-based optimum 
design of structures (consisting of a number of 
components) subjected to a '’statistical load” and 
also explicitly formulates the problem of optimi- 
zation in a framework of general cost-effectiveness 
approach. Furthermore, the present discussion on 
the statistical confidence of strength as well as load 
distribution [parts Z and 4 of this section] is more 
complete. 

6. Obviously, the present analysis is valid 
only for those cases where the quasi- static struc- 
tural analysis can reasonably well replace the 
dynamic analysis, as exemplified by the structural 
response analysis of a spacecraft to the dynamic 
pressure which builds up as a function of time in 
such a way that it will produce no significant 
dynamic effect. 

Also, in the present paper, it is assumed 
that the structural analysis accurately describes 
the stress and/or strain within the structure. The 
consideration for the error in the structural 
analysis is beyond the scope of this paper. 

VII. Conclusion 

An approach in structural optimization based 
on reliability analysis is presented with an empha- 
sis on the use of proof- load test. In particular, 
methods of optimizing the structural weight sub- 
ject to a constraint on the expected cost is 


Numerical examples, with a particular but 
reasonable expression for the expected cost, indi- 
cate that the expense of performing the proof- load 
test is always well compensated by the improve- 
ment of structural reliability due to such a test. 

In fact, under the constraint of the same expected 
cost, significant weight savings can be expected of 
a structure with proof- load-tested components, 
compared with the optimum weight of the structure 
consisting of components that are not proof- load- 
tested. The extent to which such an extra weight 
saving can be achieved depends on a parameter 
pertaining to the importance of individual com- 
ponents relative to the cost of failure. 

Also, as long as optimum levels of proof loads 
turn out to be within a central portion of the strength 
distribution, the proof-load test improves confi- 
dence of the estimated reliability value of the struc- 
ture, since if this happens, the confidence in the 
reliability estimation depends mainly on the accu- 
racy of the load prediction. An approach, in which 
such a restriction on the optimum level of proof 
load for the analysis to be meaningful can auto- 
matically be eliminated, is described. 

The question of how to deal with the statistical 
confidence of the load distribution is also discussed 
at length. 
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Appendix 

Since at optimum 6W = 0 whereas SEC'*' = 0 by 
constraint. 



n 

b.&A. - 0 

(A- 1 ) 


i- 1 



n 

E 

i=l 

' 1=1 

❖ 

aEc. 

-3X-2 6A. = 0 
dA. 1 

l 

(A- 2) 
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Hence, 


i=l 


dEC? / 0EC*\ 

l "5sf 6 'i + \ b i + x -asf J 


6 A . 


= 0 


A sufficient condition for Eqs. (A-4)[or Eqs. 
(A-6)Jto be satisfied, given Eq. (A-7) for Eq. 

(A- 8)], is 


EC* W, 

X 1 



(A-9 ) 


from which it follows that 


$ 

8EC. 

~ai7~ 


b. + X 
1 


s}c 

8EC. 

X 

dA. 

1 


0 


0 


for i = 1,2, n (A- 3) 


for i = 1,2,*** ,n (A-4) 


This statement can be verified as follows. If 
Eq. (A-9) is valid. 



and, hence, according to Eq. (A- 7) 


(A- 10) 


where X is the Lagrange multiplier. 
Equations (A-4) reduce to 


9EC* dEC* 

wj" = 


8EC 


= aW 


Therefore, 

(A- 5 ) 


and hence 



AEC* 

AEC* 



(A- 11) 


AEC? AW. 


Equations (A-6) are automatically satisfied because 
of Eqs. (A- 8), (A-9), and (A- 11). Q.E.D. 

(A-6) 
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in Eq. (1), for which Eqs. (A-3) are satisfied] 
that a small change in vi (which is proportional to 
W^) results in a significant change in EC* and 
similarly a small change of EC* results in a 
negligible change in v|. Therefore, it is reason- 
able to assume that the variation of W^/SW* due to 
that of ECj is approximately zero (the smaller the 
value of EC*, the better this approximation will 
be), i. e. , 


{w)-o ,A - 7 > 


W. + AW. W. 
. 1 1 - 1 
* • W + AW ” W 


Hence 


AW. W. 
W = “W 


Hilton, H. and Feigen, M. , "Minimum Weight 
Analysis Based on Structural Reliability," 

J. Aerosp. Sci . , Vol. 27, pp. 641-652, I960. 

^Kalaba, R. , "Design of Minimal-Weight Struc- 
tures for Given Reliability and Cost, " J. Ae rosp . 
Sci . , pp. 355-356, March 1962. 

3 

Switzky, H. , "Minimum Weight Design with 
Structural Reliability, " pp. 316-322, AIAA 5th 
Annual Structures and Materials Conference, 

New York, 1964. 

4 

Broding, W. C. , Diederich, F. W. , and Parker, 
P. S., "Structural Optimization and Design Based 
on a Reliability Design Criterion, 11 J. Spacecraft , 
Vol. 1, No. 1, pp. 56-61, January 1964. 

^Murthy, P. N. and Subramanian, G. , "Minimum 
Weight Analysis Based on Structural Reliability, " 
AIAA J . , Vol. 6, No. 10, pp. 2037-2038, 

October 1968. 

^ Moses , F. and Kinser, D. E. , "Optimum Struc- 
tural Design with Failure Probability Con- 
straints , " AIAA_J. , Vol. 5, No. 6, pp. 1152- 1158, 
June 1967. 


16 



7 

Freudenthal, A. M. , Garrelts, J. M. and 
Shinozuka, M. , "The Analysis of Structural 
Safety, 11 J. Structural Division , Proc. ASCE, 
Vol. 92, No. ST 1, pp. 267-325, February 1966. 

Q 

Shinozuka, M. and Itagaki, H. , "On the Relia- 
bility of Redundant Structures," Annals of Relia - 
bility and Maintainability , Vol. 5, pp. 605-610. 

^Zoutendijk, G., Methods of Feasible Directions , 
Elsevier Publishing Co. , Amsterdam, I960. 

10 Timoshenko, S. and Woinowsky-Krieger, S. , 
Theory of Plates and Shells , McGraw-Hill Book 
Co., Second Ed. , pT 54^, New York, 1959. 


U Wang, H. C. , and Worley, W. J. , "An 

Approach to Optimum Shape Determination for 
a Class of Thin Shells of Revolution, n Jour, of 
Applied Mechanics , pp. 524-529, Sept. .1968. 


12 Hahn, G. H. , and Shapiro, S. S. , Statistical 
Models in Engineering , John-Wiley and Sons, 
New York, p. 220, 1967. 


^Barnett, R. L. and Hermann, P. C. , "Proof 
Testing in Design with Brittle Materials," 

J. Spacecraft and Rockets , Vol. 2, No. 6, 
pp. 956-961, Nov. -Dec. , 1965. 


17 



